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Abstract 

I consider whether we can significantly improve the Cardelli et al. (19891) family of extinction 



laws using new data and techniques. There are six different aspects that need to be treated: 
The use of monochromatic quantities, the three different wavelength regimes (NIR, optical 
and UV), the sample, and the photometric calibration. Excluding the behavior in the NIR 
and UV, I discuss the other four aspects and propose a new family of extinction laws derived 
from VLT/FLAMES and HST/WFC3 data. 

1 The CCM family of extinction laws 



The Cardelli et al. (1989) or CCM family of extinction laws are the most oft-cited 
extinction laws in the astronomical literature. Their fame is undoubtedly well deserved 
since they were the first laws that accurately described the full NIR-optical-UV range with 
a parameterized form that allowed different tyes of extinction to be considered. Almost a 
quarter of a century later they are still broadly used in a wide range of applicatons. In 
this contribution we ask whether the CCM laws are due for an update using new data and 
techniques. 

The CCM extinction laws are a single-parameter family of functions that extend from 
x = 1/ X = 0.3 /im" 1 (33 333 A) to x = 10 /mi -1 (1000 A). The parameter that characterizes 
the family is Ry = Ay/E(B — V), which has typical values close to 3 but which in some 
environments can be slightly lower or significantly higher. The CCM extinction laws are 
divided into four wavelength ranges, the NIR (x = 0.3-1.1 /xm _1 ), the optical (x = 1.1-3.3 
//m _1 ), the UV (x = 3.3-8.0 ^m _1 ), and the FUV (x = 8.0-10.0 /mi -1 ). Each range uses a 
different functional form but the laws are continuous and differentiable at each boundary. The 
extinction laws in the NIR and optical were derived from multiband ground-based photometry 
and in the UV and FUV from IUE spectrophotometry. Three examples of CCM extinction 
laws are shown in Figs, [I] and [2] 
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Figure 1: CCM and new extinction laws for three values of -R5495 (2.5, 3.2, and 5.0). Ex- 
tinction is normalized to the value at 10 000 A in each case to emphasize that the extinction 
laws are the same for longer wavelengths and to better visualize the differences in the op- 
tical and NUV ranges. The approximate extent of some filters in three common systems 
(Johnson-Cousins, Strdmgren, and WFC3) is shown. 



There are several issues with the CCM laws that need to be discussed in other to study 
the need for a new family of extinction laws: 

1. Monochromatic quantities. Using Ry as the parameter for the type of ex- 
tinction is ill defined. Ry is a filter- integrated quantity (using Johnson B and V), not a 
monochromatic on^J That means that unless A(X) is constant over the extent of the filter 
(which is never the case), Ry is not only a function of the extinction law but also of the 
input spectral energy distribution and of the amount of extinction. In other words, the same 
amount and type of dust in front of two stars produces different values of Ry and doubling 
the amount of dust without changing its type also changes Ry. For the same reason, using 
Ay or E{B — V) to parameterize the amount of extinction is also ill defined. Graphical 
examples of the problem are shown in Fig. [3j 

1 A quantity will be called "monochromatic" when it is defined based only on discrete wavelengths. 



Maiz Apellaniz 

20000 10000 



X (Angstroms) 

5000 



3 



0.3 



0.2 



0.1 



< 
< 



0.0 



-0.1 



-0.2 



0.5 



K5495 = 2.5 CCM 
«549 5 = 2.5 new 



j = 3.2 CCM 
, = 3.2 new 



"5495 " 

R S4S 5 = 5.0 CCM 
^5495 = 5.0 new 
Johnson-Cousins 
Strbmgren 
WFC3 



F814W 




F555W 



1.0 



1.5 



2.0 



2.5 

: 1 IX ((irrf 



3.0 



3.5 



4.0 



Figure 2: As Fig. [I] but with a different normalization to emphasize the differences between 
extinction laws. In each case we have subtracted a linear fit A(X) / A(5495A) = 0(^5495) + 
6(i?5495)a;, with 0(^5495) and 6(^25495) calculated so that the CCM law for that -R5495 is 0.0 
at x = 1.1 and at x = 3.3, the limits for the optical range in CCM. 



So what is going on here? Actually, a CCM law with a given Ry parameter does not ac- 
tually produce extinction with Ay /{E(B — V) equal to that value of Ry (see Fig. [3]) , an effect 
that has created confusion in the past. The easiest way to fix it is to select two wavelengths 
near the center of the B and V filter and substitute the Ry (type of extinction) and E(B — V) 
(amount of extinction) parameters by monochromatic equivalents. Maiz Apellaniz| ( 2004| ) did 
this choosing 4405 A and 5495 Aas reference wavelengths for B and V, respectively. Those 
are also the choices in the current version of CHORIZOS (v 3.2, Maiz Apellaniz 2004) and 
here we follow the same convention. The reason for using those wavelengths to define the 
equivalents of Ry and E(B — V), e.g. -R5495 and £^4405 — 5495), is that for hot stars and 
low extinction [£(4405 - 5495) < 1.0], Ry « #5495 and E(B -V) » £(4405 - 5495) (Fig. [3]). 
Therefore, we will use -R5495 as the parameter for both the CCM and new laws. 
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2. The NIR. In that wavelength range (called simply infrared in the original paper), 
CCM laws use a power law with a fixe of] exponent ^(A)/^(10 00 A) = x . However, their 
analysis has been put into question by Nishiyama et al. (2009) and different authors (e.g. 



Moore et al.|2005 Fitzpatrick h Massa|2009 ) have found that the NIR extinction law changes 



from one sightline to another, with exponents as low as 1.1 and as high as 2.3. Beyond the K 



band, the extinction law first flattens and then develops complex structures (Roman-Zuniga 



et al., 2007 Nishiyama et al., 2009 Gao et al., 2009). 



3. The optical range. In order to fit the extinction law there, CCM used a seventh- 
degree polynomial. That functional form has the advantage of being able to fit the law 
through their five passbands (UBVRI) but the disadvantage of possibly introducing unde- 
sired wiggles in wavelength (Fig. [2]). Such wiggles can be detected when comparing CCM pre- 
dictions with spectrophotometry or with intermediate-band photometry such as Stromgren's. 
Prior to CCM, Whitford (1958) (see also Ardeberg & Virdefors 1982) had used spectropho- 
tometry to propose a simpler functional form of two straight lines joined near x = 2.25 /zm _1 . 
The point where the two lines are joined is called the "knee" and is the most prominent wig- 
gle in the CCM laws in the optical (Fig. [2j see also |Stebbins & Whitford"l|1943[ ). CCM were 
clearly aware of the issue and mentioned that "The Whitford law may therefore be more 
accurate for the diffuse ISM near x ~ 2.25 fj,m~ . The virtue of ours, however, is that it joins 
smoothly onto the UV extinction law from the FM Fitzpatrick & Massa 1988 sample of 
stars and that it takes into account the differences in the extinction laws of lines of sight with 
various values of Ry" . In other words, they were willing to sacrifice accuracy (in the form 
of detailed behavior for small wavelength scales) in favor of functionality (in the form of the 
addition of the Ry parameter and the continuity and differentiability for all wavelengths) . 

4. The UV and FUV. The emphasis of CCM laid on these wavelength ranges, as 
appropriate at the end of a decade where IUE revolutionized the study of extinction. Special 
importance was paid to the ability of a single-parameter family to describe extinction from 
the IR to UV, as evident in the title of the paper. However, this aspect has been challenged 
by Fitzpatrick & Massa (2007), who find that "With the exception of a few curves with 
extreme (i.e., large) values of R(V), the UV and IR portions of Galactic extinction curves 
are not correlated with each other." 

5. The CCM sample. CCM used a sample of 29 stars, which was a relatively large 
number for the time. However, more modern extinction studies use hundreds or thousands of 
stars. Also, their sample had only low or moderate extinction: only three stars have values 
of (E(B — V) greater than 1.0 and even the most reddened (HDE 229196) of those has an 
E{B — V) of 1.22. Also, the three most reddened stars have a relatively narrow range of 
-R5495 between 2.6 and 4.2. The solution for high values of -R5495 (> 5.0) depends heavily on 
a single star, Herschel 36, which has an E{B — V) of 0.89 (the other cases with -R5495 > 5.0 
have E(B — V) < 0.4) and has been recently found to be a multiple system with nearby IR 



sources (Arias et al. 2006 2010). Using stars with low extinction to determine an extinction 



law is quite dangerous since random and systematic (see below) uncertainties can introduce 
significant errors. One should also take into account that the CCM law has been applied to 



2 This is not obvious in Fig. 3 of CCM because the extinction law there is normalized with respect to V or, 
more properly, 5495 A. 
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Figure 3: R v (left) and E(B -V)- £(4405 - 5495) (right) as a function of £(4405 - 5495) 
for a CCM extinction law with -R5495 = 3.2 and three main-sequence stars with different T e g. 
Ry ~ -R5495 and E(B — V) ~ £(4405 — 5495) only for hot stars with low extinctions. 



objects with £(4405 — 5495) significantly larger than 1.0: in those cases one is extrapolating 
a law beyond the range for which it was derived, thus amplifying any potential errors in the 
original work. 

6. Photometric calibration. The final issue is that of the photometric calibration 
of the filters used in the NIR and optical ranges. By calibration we mean both the determi- 
nation of the zero points and the shape of the sensitivity curves. The reader is referred to 



Mafz Apellaniz ( 2005 2006 , 2007 ) and references therein for developments in this field after 



CCM. In particular, it should be noted that old observations in the Johnson U band may be 
suspect because of the effect of the atmosphere in the sensitivity curve and because the curve 
straddles the Balmer jump, thus being affected by errors in the spectral type. 



2 The new extinction laws 

Based on the issues discussed in the previous section, we decided to attempt the cal- 
culation of a new family of extinction laws. Ideally, to complete such a task one would use 
high-quality spectrophotometry from the NIR to the UV of a diverse collection of sources 
in different environments and with different degrees of extinction. Since such dataset is not 
currently available, we will start by dealing only with some of the issues discussed in the 
previous section. More specifically, we will ignore extinction in the UV (except for the region 
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closest to the optical) and for the NIR we will simply use the CCM laws. In other words, 
we will concentrate our efforts in the optical region. Ignoring the UV will not matter to a 
non-specialist interested only in eliminating the extinction from his/her optical data. Ignor- 
ing the NIR may matter if the exponent there is significantly different from the CCM one 
but only if extinction is very large and even then it may only apply to the total extinction 
correction, not to e.g. the determination of T e g from the photometry. 

To derive the new family of extinction laws we are using two datasets: 



The VLT/FLAMES Tarantula Survey (Evans et al. 2011| ), from which we are using 
~ 200 accurate spectral types of stars (mostly of O type) in 30 Doradus as well as NIR 
photometry. 



• The WFC3 Early Release Science HST/WFC3 images of the central region of 30 
Doradus (De Marchi et al. 2011; Sabbi et al.[ |2012 ), from which photometry in six 
broadband filters (F336W, F438W, F555W, F814W, F110W, and F160W, which are 
UBVIJH equivalents) was extracted. 



Our analysis is in a quite advanced stage and in the near future we will submit for 
publication our results. A preliminary view with three examples of the new extinction laws is 
shown in Figs. [T] and [2] We have verified that the new laws provide significantly lower residuals 
when assuming effective temperatures derived from the spectral types and better fits when 
attempting to derive the effective temperatures from the photometry. We have also applied 
the new extinction laws to Galactic stars with £"(4405 — 5495) between 1.5 and 2.0 and have 



been able to significantly lower the fits residuals when compared to CCM or to Fitzpatrick 



(1999), thus extending the applicability of the new family of optical/NIR extinction laws to 
the Milky Waj0 
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3 Note that the extensive literature on the differences between the SMC, LMC, and MW extinction laws 
refers to the UV regime, not to the optical/NIR. 
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